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We study a two-dimensional spin-orbit-coupled dipolar Bose-Einstein condensate with repulsive 
contact interactions by both the variational method and the imaginary time evolution of the Gross- 
Pitaevskii equation. The dipoles are completely polarized along one direction in the 2D plane so 
as to provide an effective attractive dipole-dipole interaction. We find two types of solitons as the 
ground states arising from such attractive interactions: a plane wave soliton with a spatially varying 
phase and a stripe soliton with a spatially oscillating density for each component. Both types of 
solitons possess smaller size and higher anisotropy than the soliton without spin-orbit coupling. 
Finally, we discuss the properties of moving solitons, which are nontrivial because of the violation 
of Galilean invariance. 
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I. INTRODUCTION 

Ever since the first achievement of Bose-Einstein con¬ 
densates (BECs) in ultracold atomic gases [1], matter 
wave solitons have been the central focus of many exper¬ 
imentalists and theorists [2] . Solitons are the result of the 
interplay between nonlinearity and dispersion and keep 
their shape while traveling. In BECs, nonlinearity origi¬ 
nates from collisional interactions between atoms, which 
can be readily tuned via Feshbach resonances [3]. In gen¬ 
eral, there are two types of solitons: a bright soliton with 
a density bump for attractive interactions and a dark 
soliton with a density notch and a phase jump for re¬ 
pulsive interactions. Both bright and dark solitons have 
been experimentally observed in cold atoms with contact 
interactions [4-12]. However, for such contact attractive 
interactions, bright solitons can only exist in one dimen¬ 
sion (ID), but not in two dimensions (2D) where the 
states either collapse or expand [13]. 

Different from the local nonlinearity resulting from 
contact interactions, the non-local nonlinearity can sta¬ 
bilize a 2D bright soliton [14, 15], in particular, the 
nonlinearity introduced by the dipole-dipole interaction. 
This interaction is long ranged and anisotropic with the 
strength and sign (i.e. repulsive or attractive) depend¬ 
ing on the dipole orientation. When an external rotating 
magnetic field is applied to reverse the sign of the inter¬ 
action [16], or the dipoles are completely polarized in a 
2D plane [17], the dipolar interaction can become attrac¬ 
tive and 2D bright solitons can be, therefore, stabilized 
under appropriate conditions. It is essential to note that 
although the relevant interaction in common experiments 
with cold atomic gases is contact, increasing interest has 
been focused on the atoms with large magnetic moments 
that possess dipole-dipole interactions [15, 18, 19]. In 
fact, the Bose-Einstein condensation of several dipolar 
atoms such as Chromium [20-22], Dysprosium [23], and 
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Erbium [24] , as well as the degeneracy of a dipolar Fermi 
gas [25, 26] have been observed in experiments. 

Recently, the spin-orbit coupling between two hyper- 
fine states of cold atoms has been experimentally en¬ 
gineered [27-32]. And this achievement has ignited 
tremendous interest in this field because of the dramatic 
change in the single particle dispersion (induced by spin- 
orbit coupling) which in conjunction with the interaction 
leads to many exotic superfiuids [33-43] (also see [44-51] 
for review). Such change in dispersion also results in 
exotic solitons even when the interaction is contact, in¬ 
cluding bright solitons [52-58], dark solitons [59, 60], and 
gap solitons [61-63] for BECs, as well as dark solitons for 
Fermi superfiuids [64, 65]. These solitons exhibit unique 
features that are absent without spin-orbit coupling, for 
instance, the plane wave profile with a spatially varying 
phase and the stripe profile with a spatially oscillating 
density for BECs, as well as the presence of Majorana 
fermions inside a soliton for Fermi superfiuids. Also, the 
violation of Galilean invariance [53, 66, 67] by spin-orbit 
coupling dictates that the structure of solitons changes 
with their velocities. 

On the other hand, spin-orbit-coupled BECs with 
dipole-dipole interactions [68-71] have also been ex¬ 
plored, and intriguing quasicrystals [72] as well as meron 
states [73] have been found. However, whether a soliton 
can exist in such BECs in 2D with long ranged dipole- 
dipole interactions and spin-orbit dispersion has not yet 
been investigated. 

In this paper, we examine the existence and properties 
of a bright soliton in a two species spin-orbit-coupled 
dipolar BEC in 2D with repulsive contact interactions 
via both the variational method and the imaginary time 
evolution of the Gross-Pitaevskii equation (CPE). The 
dipoles are completely oriented along the y direction in 
the 2D plane in order to provide an effective attractive 
dipole-dipole interaction. Thanks to such attractive in¬ 
teractions, we find two types of solitons: a plane wave 
soliton (when the repulsive intraspecies contact interac¬ 
tion is larger than the repulsive interspecies one) and a 
stripe soliton (when the interspecies one is larger). These 
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2D solitons as the ground states cannot exist for a system 
with pure attractive contact interactions and spin-orbit 
coupling. Such solitons are highly anisotropic and their 
size is also reduced by spin-orbit coupling. Finally, we 
study the moving solitons, which are nontrivial because 
of the lack of Galilean invariance. The size of a soli- 
ton hrst increases and then decreases with the rise of the 
velocity and this change is anisotropic. The moving soli- 
ton also tends to be plane wave even when its stationary 
counterpart has the stripe structure. 

The paper is organized as follows. In Sec. II, we intro¬ 
duce the energy functional and the time-dependent GPE, 
which are used to describe a spin-orbit-coupled dipolar 
BEG. In Sec. Ill, we calculate the bright soliton by per¬ 
forming the minimization of the energy of the variational 
ansatz wave functions and an imaginary time evolution 
of the GPE. The properties of such soliton are also ex¬ 
plored by the former method. Then, we study the non¬ 
trivial moving solitons in Sec. IV. Finally, we conclude 
in Sec. V. 


II. MODEL 


states) of a large spin atom (e.g. dysprosium) is consid¬ 
ered. We also assume that the dipoles are all oriented 
along the y direction, thus 

Edd= J ( 3 ) 


where the Fourier transform of the total density is pk = 
/ -|- |d>|.p) and C/(k) is given by 


t/d(k/.) 


. — Snlzkle^ ^*/^erfc(fcG/A/2) 

V^-b-^^^ 

k 


(4) 


with erfc being the complementary error function. Here, 
gd = characterizes the strength of the 

dipole-dipole interaction where p is the magnetic dipolar 
moment and po is the permeability of the free space. 

The dynamical behavior of a BEG can be described by 
the time-dependent GPE 




Hs'i’ir) + pGT(r) -b gdUd{r)'if{r), 


(5) 


where the contact interaction matrix is 


We consider a Rashba-type spin-orbit-coupled BEG 
and write its single particle Hamiltonian as 


G = 


0 



( 6 ) 


P‘ , 2 2 1 


-ffs = ^ + -b A(p x cr) • e,,, ( 1 ) 


where p = —ifiV is the momentum operator, m is the 
atom mass, A is the spin-orbit coupling strength, and cr 
are Pauli matrices. uj± (uJz) is the trap frequency in the 
(x,y) plane (along the z direction). Here, we assume 
that hcJz is much larger than huj± and the mean-field 
interaction so that the atoms are frozen to the ground 
state in the z direction. Given that a soliton is studied, 
we thus set w_l = 0. 

When the s-wave contact and dipole-dipole interac¬ 
tions are involved, the energy functional of a 2D conden¬ 
sate can be written as 





+5i2|4'-|-p|d>j_P] + Edd, 


( 2 ) 


where the condensate wave function d'(r) = 

[^'■^(r), d' 4 ^(r)]^ with two pseudo-spin compo¬ 
nents 'I'.|.^)(r) , g = Airh^a /and 
gi 2 = EKk^ayi !are the intraspecies and 
interspecies contact interaction strength respec¬ 
tively with the intraspecies and interspecies s-wave 
scattering length being a and a \2 and the charac¬ 
teristic length along z being = \Jh,l (muiz)- Here, 
Hg = + dy)/{2m) — ihX{dxCry — dycr^) is the 2D 

single particle Hamiltonian, and the third dimension 
has been integrated out. For dipole-dipole interactions, 
we only consider the density-density interaction which 
is dominant when a two subspace (i.e. two pseudo-spin 


and the dipolar interaction potential is 

1 


Ud{Y) = 


(27r)2 


dke*‘^-Xk)C/d(kG). 


(7) 


For numerical simulation, we choose tvjJz, Iz, and l/wz 
as the units of energy, length, and time, respectively, and 
the dimensionless energy per atom hence reads 


£= dr 


$(r)t7^,$(r) + i7(|$^|4 + |ci>^|4) 


+7i2i‘f’tn4^f n + 


Id 


2(27r)2 


dknkn_kC/d(k), (8) 


where 'Hs = -{dl + i9y)/2 - ia{dxcry - dyax), a = 
X/{uJxlz), 7 = 2y/^Noa/lz, 712 = 2 y/^NQai 2 /lz with 
the total particle number iVo, , 7 ^ = 2NQad/lz with 
ad = mkok?/{12-Kh?), andrik = / (ire“*'^''’(|$i-P + |4>4,P)- 
The wave function is normalized to 1 (i.e. /(ir(|<I>.|-p -b 
|4>iP) = l). 

The dimensionless time-dependent GPE reads 

i^^=H.$(r)+7a$(r) 


dt 


+ ■ 


Id 


where 


Q = 


(27r)2 

0 


J dkE^ '^nik)Udik)<i>{r), (9) 




( 10 ) 
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FIG. 1. (Color online) Profiles of the density of t O') in (a-)((b)), the total density + nj. in (c), the phase 

of spin 'I'O') in (d)((e)) for a plane wave soliton (the first two panels) with 712 = 6 and a stripe soliton (the last two panels) 
with 712 = 10. The solitons in the first and third panels are obtained by the variational method, while the solitons in the 
second and forth panels are calculated by the imaginary time evolution of the GP Eq. (9). The dashed white line labels the 
a; = 0 line. Here, 7 = 8 , 7^/7 = 0.67, and a = 2. 


III. STATIONARY BRIGHT SOLITONS 

To shed light on the structure of a soliton, we start 
from the homogeneous noninteracting single particle sce¬ 
nario and write its momentum space dispersion as 

£^(k) = Y±afc, (11) 

with two branches labeled by the helicity ±. Clearly, the 
ground state is degenerate with the energy being — 0^/2 
when the momenta lie in the k = \a\ ring. This is differ¬ 


ent from the case without spin-orbit coupling where the 
ground state only occurs at fc = 0. In this single particle 
case, any superposition of the states in the ring is also its 
ground state. Yet, this is not the case when the repulsive 
contact interaction is involved. The ground state either 
possesses a single momentum (i.e. plane wave phase) 
when 712/7 < 1 or two opposite momenta (i.e. stripe 
phase) when 712/7 > 1 [34]. When the dipolar inter¬ 
action is turned on, one may expect that this effective 
long ranged attractive interaction along with contact re¬ 
pulsive interaction could support two types of solitons: 
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FIG. 2 . (Color online) Plot of in (a) and ay in (b) as 
a function of for the plane wave solitons (dotted blue 

line), stripe solitons (dashed green line), and traditional soli¬ 
tons (solid red line) without spin-orbit coupling. The aspect 
ratio ^/ayjax, of a soliton is displayed in the inset of (b). 
Variational parameters with respect 7d/7 are plotted in (c) 
associated with xo (dash-dot blue line) and Jp (dashed blue 
line) for the plane wave soliton, and Jx (solid green) and Jy 
(dotted green) for the stripe soliton. In (d), the total energy 
of the variational ansatz wave function and the wave function 
numerically obtained by the imaginary time evolution of the 
GPE for both plane wave and stripe solitons is shown. The 
solid green line (stripe soliton) and dashed blue line (plane 
wave soliton) correspond to the variational results, while the 
green circles and blue squares correspond to the GP results. 
Here, a = 2 , 7 = 8, 712 = 6 (712 = 10 ) for the plane wave 
(stripe) soliton. 

plane wave and stripe solitons. 

To examine whether a soliton can indeed exist in 
the spin-orbit-coupled dipolar BECs, we first consider a 
plane wave soliton variational ansatz 

‘I’P = exp(-zJpy), (12) 

where 

V 27r 

Here Jp is the wave vector of the plane wave soliton, = 
1/v^^ with V = x,y IS the size of the soliton, and cco is 
the separation distance between two components. When 
xq = 0, this state is an eigenstate oipyUx multiplied by a 
Gaussian profile <f>o(0), and Jp = a yields the minimum 
energy. In fact, Xq is usually nonzero because of a force 
acting on the BEG by spin-orbit coupling F = a^(p x 


ez)o'z [38, 74], which is opposite along the x direction 
when p (here Jp) is along the y direction. 

In writing down the ansatz (12), we have assumed that 
the wave vector Jp is in the y direction. The prerequisite 
of this assumption is that the rotation symmetry [38, 39] 
about the z axis has been broken by the dipole-dipole 
interaction. Indeed, without the dipole-dipole interac¬ 
tion, this state with Jp along y is not special and other 
states with Jp along other directions are degenerate with 
it. For example, the state with Jp along y has the same 
energy as a state with Jp along x. Yet, with the spe¬ 
cific dipole-dipole interaction arising from the dipoles en¬ 
tirely oriented along y, the symmetry is broken and the 
ground state should be elongated along y (a^ > ay) so as 
to provide an effective attractive interaction because of 
the head-to-tail configuration of polarized dipoles. This 
elongated configuration allows the existence of a 2D soli¬ 
ton [17] and also requires the wave vector to be along 
y [75]. 

Although the wave vector Jp of the ground state is 
along y, there are still two options: negative and positive 
directions in terms of the time-reversal symmetry "T (i.e. 
—icjylC with the complex conjugate operator /C). Specif¬ 
ically, the state <I>P 2 =T^p is degenerate with $p. In 
the absence of interactions, all superposition states of $p 
and 4 >p2, 

‘hps = I cos01$p-b 1 sin0je*‘^$P2, (14) 

are degenerate. This degeneracy may be broken by the 
interaction so that the ground state is either $p or <I>P2, 
or a certain superposition state of them. But this de¬ 
generacy breaking should not happen at 712/7 = 1 since 
the interaction energy only depends on the total density 
which is independent of 6 and p. This gives us an in¬ 
tuitive understanding that 712/7 = 1 may separate the 
plane wave soliton (j cos 9\ =0 or 1) and the stripe soliton 
(1 cos^l = 1 sin0j), similar to the homogenous spin-orbit- 
coupled BEG [34] without dipole-dipole interactions. For 
the stripe soliton, we note that (p = 0,7r corresponds 
to the ground state as the energy contributed by p is 
- 7 i 2 v/^e-(“-“«/ 2 + 2 J=/aD cos(2^)/(167r) [76]. 

To evaluate the variational parameters Qx, ay, Xq, Jp, 
and 9, we minimize the energy £ after substituting <I>ps 
in Eq. (14) to Eq. (8). Indeed, the calculated variational 
solutions reveal that there are two types of soliton so¬ 
lutions: plane wave solitons when 712/7 < 1 and stripe 
solitons when 712/7 > 1. We present the density and 
phase profiles of a typical plane wave soliton (we choose 
9 = tf = Q) in the first panel of Fig. 1, where the stripe 
structure of the phase of both two components reveals 
the plane wave feature. The soliton is highly elongated 
along the y direction and the centers of two components 
are spatially separated along the x direction because of 
nonzero xq. To confirm that this variational solution can 
qualitatively characterize the ground state of the system, 
we numerically compute the ground state by an imagi¬ 
nary time evolution of the GP Eq. (9). This exact nu¬ 
merical solution also concludes that 712/7 < 1 yields the 
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plane wave soliton while 712/7 > 1 the stripe soliton. 
In the second panel of Fig. 1, we also plot the corre¬ 
sponding density and phase profiles of the GP obtained 
plane wave soliton. The variational ansatz is in qualita¬ 
tive agreement with it given the separated centers and 
the plane wave varying phase that both states possess. 
Yet, the shape of the soliton obtained by the imaginary 
time evolution deviates slightly from the Gaussian and 
the size is also slightly smaller. 

When 6 = 7r/4 and (/? = 0, ^ps is a stripe state 
with a density oscillation along the y direction for each 
component. And there is no stripe for the total den¬ 
sity. Along the x direction, two components are not spa¬ 
tially separated, and the phase for the spin f reverses 
suddenly across x = 0. Following these properties by re¬ 
placing ($o(a;o/2) -I- ^oi-xo/2))/V2 with cos{JxX)^o{0) 
and (4>o(xo/2) — $o(—a;o/2))/-\/2 with sin( J 2 ,x)<I>o( 0 ) in 
Eq. (14), we obtain another better variational ansatz for 
the stripe soliton 


$s=r$o(0), (15) 


where 

T = ( cos(Ja;x) - i sm{Jyy) sin(J,;,a;)\ . 

\cos{Jyy) sin(Ja;a;) -I- f sin(Jy?/) cos{Jxx)J ’ ^ ' 

with the variational parameters Jx and Jy. The period 
of the stripe along the y direction is tt/J y. Interestingly, 
this stripe soliton corresponds to four points (± J^, ± Jy) 
in momentum space instead of traditional two points [53] 
when Ja; = 0, if we do not consider the Gaussian profile 
$ 0 . 

We calculate the variational parameters of stripe soli- 
tons by performing the minimization of the energy £ in 
Eq. ( 8 ) where $ is replaced with $ 5 . The density and 
phase profiles of a typical stripe soliton calculated by this 
method is displayed in the third panel of Fig. 1. Evi¬ 
dently, the density of each component exhibits the stripe 
structure while the total density does not. The phase 
of spin t along the y direction varies like a plane wave, 
but reverses across a; = 0 due to the presence of sin( Ja;x) 
in the imaginary part. The phase of spin J, exhibits the 
phase rotation like vortices around a; = 0 and y = mr / Jy 
with integer n; around these points, the wave function 
$ 54 , is proportional to (—1)"( J^a;-I- 1 ( Jy?/ — nTr)) and the 
corresponding density of spin is extremely low. More¬ 
over, in the last panel of Fig. 1, we plot the density and 
phase profiles of the corresponding stripe soliton obtained 
by the imaginary time evolution of the GPE; comparing 
this figure with the third panel of Fig. 1 implies that the 
stripe variational ansatz is qualitatively consistent with 
the GP results. 

To study the properties of a soliton with respect to 
dipole-dipole interactions 7 ^, we evaluate the variational 
parameters of both the plane wave and stripe solitons by 
the variational method and plot them in Fig. 2 as 7^/7 
varies. Clearly, with increasing 7 d/ 7 , a-x and Oy increase 
monotonously because of the enhanced effective attrac¬ 
tive interaction, indicating that the size lx and ly of the 



FIG. 3 . (Color online) Plot of xo (dash-dot blue line) and Jp 
(dashed blue line) for the plane wave variational ansatz, and 
Jx (solid green line) and Jy (dotted green line) for the stripe 
variational ansatz with respect to by the minimization of 
the energy and in (a), (b) shows the minimum 

energy of Af(dotted blue line) and (solid green line) 

as a function of a^. Here, a = 2 . 


soliton decrease monotonously. We note that as 7^/7 
increases further, the soliton can collapse so that both 
ttx and Uy diverge. For the plane wave soliton, ax and 
fly are slightly larger than the stripe soliton because of 
the smaller contact interaction of the former. Moreover, 
compared with the soliton without spin-orbit coupling 
(red line in Fig. 2(a) and (b)), Qx and Oy for both the 
plane wave and stripe solitons are much larger, implying 
that the size of solitons can be reduced by spin-orbit cou¬ 
pling. Also, these solitons are highly anisotropic with the 
much smaller aspect ratio y'ayjox as shown in the inset 
of Fig. 2(b). To elucidate the reason, we explicitly write 
that single particle energy of the plane wave variational 
ansatz in Eq.(12) which results from the presence of xq 
and Jp 

= ^Jp-(^(^~~^{Jp + \axXo). (17) 

The minimization of with respect to xq and Jp for 
fixed ax yields 

—Jp + ^ / Jp -l- 2ax 

*0 =- (18) 

Jp = (19) 

For ax = 0, the energy is independent of xq and Jp = a, 

while for Oj, 7 ^ 0 , both Xq and Jp decrease slightly with 
increasing ax as shown in Fig. 3(a) with the asymptotic 
Xq = 1/a and Jp = a as ax goes zero. The energy 
is also a monotonously decreasing function of ax- And 

this energy decline combined with the reduced dipole- 
dipole interaction energy competes with the rise of the 
kinetic energy (when xq = Jp = 0) and contact inter¬ 
action energy, leading to an increased ax and Oy com¬ 
pared with the soliton without spin-orbit coupling. This 
is also consistent with Fig. 2(c), showing that with in¬ 
creasing the dipole-dipole interaction, ax increases and 














6 



a a 


FIG. 4 . (Color online) Plot of ax in (a) and ay in (b) as a func¬ 
tion of the spin-orbit coupling strength a for the plane wave 
soliton (dotted blue line) and stripe soliton (dashed green 
line). The aspect ratio ^fayjax of a soliton is plotted in the 
inset of (a), (c) illustrates the change of xo (dash-dot blue 
line), Jp (dashed blue line) for the plane wave soliton and 
Jx (solid green line) and Jy (dotted green line) for the stripe 
soliton with respect to a. In (d), the total energy plus o? 12 
is plotted as a function of q. The solid green (for a stripe 
soliton) and dashed blue (for a plane wave soliton) are calcu¬ 
lated by the variational method, while the green circles (for a 
stripe soliton) and blue squares (for a plane wave soliton) are 
numerically obtained by the imaginary time evolution of the 
GPE. Here, 7 = 8 , 7^/7 = 0.67, and 712 = 6 (712 = 10) for 
the plane wave (stripe) soliton. 

both Xq and Jp, therefore, decrease so as to lower Ef^. 
It is important to note that although is not a func¬ 
tion of Qy, other energy such as the kinetic energy (when 
Xo = Jp = Q), the contact and dipolar interaction energy 
depends on it. 

For the stripe soliton, the single particle energy due to 
the presence of Jx and Jy is 

= \{Jl + Jl) - a{Jx + Jye-^"!’^-). (20) 

Similar to the plane wave case, this energy is indepen¬ 
dent of ay. For fixed ax, the minimization of this energy 
yields both Jx and Jy as a function of ax as shown in 
Fig. 3(a). When ax moves towards zero, the solution ap¬ 
proaches {Jx = a, Jy = 0) or {Jx = 0, Jy = a); when it 
moves away from zero, there is only one solution where 
Jy decreases from a while Jx increases from zero with 
the rise of Qx- Also, the energy decreases as ax 

increases. Analogous to the plane wave soliton, the to¬ 
tal energy decrease resulted from spin-orbit coupling and 
dipole-dipole interactions as ax and ay increase from the 


value without spin-orbit coupling exceeds the energy gain 
of the kinetic (when J^, = 0 and Jy = 0 ) and contact in¬ 
teraction; this leads to the increased ax and ay compared 
with the soliton without spin-orbit coupling. This picture 
is also consistent with Fig. 2(c) where Jx increases while 
Jy decreases with respect to 7 ^/ 7 . 

To explicitly demonstrate the effect of the spin-orbit 
coupling on the properties of a soliton, we plot the varia¬ 
tional parameters as a function of the spin-orbit coupling 
strength a for both the plane wave and stripe solitons 
in Fig. 4. Consistent with the aforementioned feature 
that spin-orbit coupling can reduce the size of the soli¬ 
ton, both Fig. 4(a) and Fig. 4(b) display a monotonous 
increasing behavior of ax and ay as a function of a. Also, 
the aspect ratio yjayjax is decreased by spin-orbit cou¬ 
pling. Similar to Fig. 2(a) and Fig. 2(b), ax and ay for 
the plane wave soliton are slightly larger than the stripe 
soliton in that the former has a smaller contact interac¬ 
tion. For the plane wave soliton, Jp (determined mainly 
by the spin-orbit coupling strength) increases with re¬ 
spect to a while Xq decreases; for the stripe soliton, both 
Jx and Jy increase. 

In Fig. 2(d) and Fig. 4(d), for both plane wave and 
stripe solitons, we compare their energy obtained by the 
variational procedure with the one obtained by the imag¬ 
inary time evolution of the GPE. Both figures show that 
the energy calculated by the imaginary time evolution is 
lower as expected. Yet, the difference between these two 
energy is not large (no more than 10 %), suggesting that 
the variational ansatz can qualitatively characterize the 
solitons. We note that in Fig. 4(d), the energy is shifted 
by 0^/2 in order to clearly present the different results 
of the two methods, which could be smeared by the large 
value of a^/ 2 . 

IV. MOVING BRIGHT SOLITONS 

Generally, the wave function of a moving soliton with 
the velocity v can be simply written as exp(iv • r)$s(r — 
vt) where <I>s is the wave function of a stationary soliton. 
But this is only valid for a system respecting Galilean 
transform invariance. In fact, Galilean invariance is bro¬ 
ken in a spin-orbit-coupled BEG [ 66 ], and this violation 
dictates that the shape of a soliton depends on its velocity 
strength [53]. Here, for a soliton in a spin-orbit-coupled 
dipolar BEG in 2D, we assume that a moving soliton can 
be written as 

$M(r, t) = 4>„(r - vt, t) exp(iv • r - ( 21 ) 

where is a localized function. Plugging 4>M(r,t) into 
Eq. (9) yields 

0 ^ ( r'l 

»^^^=l^«(v)<i>„(r)+7e$„(r) 

+7d J dke*‘^-^n(k)C/d(k)$,(r), (22) 
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FIG. 5. (Color online) Imbalance Imb and width lx of spin 
t of the solitons with respect to the velocity along the y (x) 
direction in (a) ((b)). The insets display the enlarged figure 
in a small velocity region. Density and phase profiles of four 
typical moving solitons for spin 7 corresponding to the dif¬ 
ferent velocities in (a) and (b) are plotted in (c-f) where the 
horizontal and vertical coordinates are x and y respectively. 
Here, a = 2, 7 = 8 , 712 = 10, and 7 d /7 = 1, corresponding 
to a stripe soliton when stationary. 


where 7^s(v) = 'Hs + a(v x cr) • e^. Compared to 
Eq. (9), this dynamical equation has an additional term 
q;(v X cr) • 02, acting as a Zeeman field; this additional 
term implies the violation of Galilean invariance. This 
violation means that it is no longer a trivial task to find 
a moving bright soliton for a BEC with spin-orbit cou¬ 
pling; we need to perform an imaginary time evolution 
of the Eq. (22), but not Eq. (9). Furthermore, such a 2D 
moving soliton should be different for different velocity 
directions even if their amplitude is the same, in contrast 
to a ID soliton which can only move in one direction. 

To examine how the shape of a soliton changes with 
respect to the velocities along x and y directions, we plot 
the imbalance Imb and the width of a soliton of spin f 
as a function of the velocities Vx and Vy in Fig. 5. Here, 
the imbalance for spin f is defined as 


|$t(0)p-|cl>^(^/(2a))p 

|$t(0)P + |cI>t(V(2a))P’ 


(23) 


which characterizes a stripe soliton (as shown in Fig. 5(c) 
and Fig. 5(d)) when it approaches one and a plane wave 
soliton (as shown in Fig. 5(e) and Fig. 5(f)) when it ap¬ 
proaches zero. Fig. 5(a) and Fig. 5(b) demonstrate that 
Imb suffers a sharp decline from one to near zero as Vx 
and Vy increase, indicating that a moving soliton tends to 
be a plane wave state. The reason is the broken rotation 
symmetry of the single particle Hamiltonian by the veloc¬ 
ity induced Zeeman field, giving rise to a ground state of 


the single particle system lying at one momentum point 
located along the x {y) direction when the velocity is 
along that direction. This also explains why the phase of 
a moving plane wave soliton with the velocity along the 
X {y) direction varies along that direction. 

Furthermore, Fig. 5(a) demonstrates that the width of 
the soliton gradually grows when the velocity along the 
y direction is enlarged. To explain the growth, we con¬ 
sider the plane wave ansatz in Eq. (12) which yields an 
additional term for the single particle en¬ 

ergy when a soliton moves; this energy decrease enlarges 
exponentially with the decline of (i.e. increase of the 
width), leading to an expanded soliton with the rise of 
the velocity. However, this is not a monotonous behavior 
and the soliton begins shrinking when the velocity goes 
larger, due to the enlarged Jp by the velocity induced 
Zeeman field, similar to increasing spin-orbit coupling. 
On the other hand, when the velocity is along the x di¬ 
rection, the width of the soliton gains a sudden rise as 
the velocity varies, as shown Fig. 5(b). This corresponds 
to a change from a stripe soliton with the wave vector 
along the y direction to a plane wave soliton with the 
wave vector along the x direction. For the stationary 
solitons, the soliton with the wave vector mainly along 
the y direction has lower energy than the one with the 
wave vector mainly along the x direction as the dipoles 
are completely oriented along y. But the Zeeman field 
induced by the presence of a velocity along the x direc¬ 
tion gives rise to the single particle ground state that 
possesses the wave vector along x. The two states with 
the wave vector along these two directions compete and 
change from the former to the latter (i.e. first order phase 
transition). For the decrease of the width when Vx goes 
even larger, the reason is the same as the case for Vy. 
When a stationary soliton is plane wave, the moving be¬ 
havior is similar except that the moving soliton is always 
the plane wave soliton. 


V. CONCLUSION 

We have studied the bright solitons as the ground 
states in a spin-orbit-coupled dipolar BEC in 2D with 
dipoles completely polarized along one direction in the 
2D plane. It is important to note that the solitons are 
the ground states in 2D, but they are the metastable 
states in quasi-2D where the true ground state would 
collapse and there is an energy barrier between the soli¬ 
ton state and this ground state. Two types of solitons 
have been found: a plane wave soliton and a stripe soli¬ 
ton. The former has the plane wave phase variation and 
its two components are slightly spatially separated; while 
for the latter, the density of each component is spatially 
oscillating and the variational ansatz suggests that four 
points in momentum space are involved. Both plane wave 
and stripe solitons are highly anisotropic and their size 
is decreased by spin-orbit coupling. These solitons can¬ 
not exist as the ground states in a 2D system with pure 































attractive contact interactions and spin-orbit coupling. 
Moreover, the shape of these solitons changes with their 
velocities due to the absence of Galilean invariance, and 
this change is anisotropic. 

The 2D bright soliton, albeit mainly plane wave soli- 
ton, can also exist when equal Rashba and Dresselhaus 
spin-orbit coupling is considered. In experiments, this 
type of spin-orbit coupling has been engineered by cou¬ 
pling two hyperfine states of atoms through two counter- 
propagating Raman laser beams [27-32] and such setup 
could be employed to realize this spin-orbit coupling in 
Dysprosium [69] with large dipole-dipole interactions. 
Also, the large magnetic moment in Dysprosium atoms 


may permit the realization of Rashba spin-orbit cou¬ 
pling [77]. 
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